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Abstract 

We explore the relation between resummation and explicit multi-loop calculations 
for QCD hard-scattering amplitudes. We describe how the factorization properties 
of amplitudes lead to the exponentiation of double and single poles at each order 
of perturbation theory. For these amplitudes, previously-observed relations between 
single and double poles in different 2 — > 2 processes can now be interpreted in 
terms of universal functions associated with external partons and process-dependent 
anomalous dimensions that describe coherent soft radiation. Catani's proposal for 
multiple poles in dimensionally-continued amplitudes emerges naturally. 
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1 Introduction 

The past few years have seen a breakthrough in the calculation of two-loop 
matrix elements and amplitudes in QCD [1,2,3,4]. At the same time signifi- 
cant progress has occured in the resummation of logarithmic corrections to all 
orders in perturbation theory [5,6,7,8]. In this paper, we hope to contribute 
toward linking advances in fixed-order QCD amplitudes and in resummation. 
We attempt to clarify how the factorization properties of QCD hard scattering 
amplitudes influence calculations at fixed order. 
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We will see how the single pole structure found at the level of two loops 
in matrix elements can be understood in terms of properties of a hitherto 
uncalculated two-loop anomalous dimension for soft radiation. The knowledge 
of this two- loop soft function is an important part of resummations at the level 
of next-to-next-to-leading logarithms. In addition, we will relate Catani's very 
successful proposal [9] for the e-pole structure of dimensionally regularized 
{D = 4 — 2e) amplitudes at one and two loops to known factorization and 
resummation properties of QCD amplitudes. These results will be illustrated 
for amplitudes involving quarks and antiquarks. 



2 The factorized hard scattering ampHtude 

Consider the partonic process 

f : fA{iA, ta) + fniiB, tb) ^ r,) + Mp2, r,) , (1) 

which involves four partons with flavors /j, momenta {ii,Pi} and color rj. 
In perturbation theory, the hard scattering amplitude M.^^'^ associated with 
this process can be written as a matrix in the space of the color indices of 
its external partons n = {^a, ^s; ri, r2}. It is convenient to express these 
amplitudes in a basis of color tensors, (c/)^.., so that [10] 



K.} (i^^>' ^' «^(/^')' = {{Pjh ^, 0^s{n'),e^ {cL){r,y ■ (2) 

Here i,j — A,B,1,2 and is the overall hard-scattering scale, which we 
may choose as s for 2 — > 2 scattering. The vectors pj = pj/ \f s/2, specify the 
directions of the incoming and outgoing lines and hence the full kinematics 
of the process ^ . After renormalization, these amplitudes retain poles in e, 
due to infrared and coUinear singularities. The singular behavior of the 2 — > 
2 amplitude is entirely due to internal lines that carry momenta that are 
either soft or collinear to one of the external lines. For gauge-invariant sets 
of diagrams, these singularities are entirely logarithmic in four dimensions, 
corresponding to poles at e = 0. 

The amplitude TVf'^l for partonic process f can be factorized into products 
of functions that organize the contributions of momentum regions relevant 
to e poles in the scattering amplitude [11,12]. These are: i) process-dependent 
functions h^P that describe the short-distance dynamics of the hard scattering, 

^ For more than two particles in the final state, we must include information on 
the relative energies of the outgoing Unes. 
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one function for each of the elements of the basis of color exchange; ii) a 
matrix of functions, Sj^j that describes the coherent soft radiation arising 
from the overall color flow and iii) a function j'^, dependent only on the 
list of external partons, and otherwise independent of the color flow, which 
describes the perturbative evolution of the incoming and outgoing partons 
for flavors fi — q, q, g. In these terms, the factorized form of the amplitude 
associated with the process f of Eq. (1) is [11,10] 



X hf[p,,^,a,{^^'')^ , (3) 

where we collect the various virtual "jet" factors associated with external 
partons in 



J i=A,BX2 ) 



(4) 



The products are over the two incoming lines, and over the two outgoing lines 
(or more, in the cases where we generahze the 2 — > 2 process of Eq. (1)). 
Notice that in Eq. (3) only the soft and hard-scattering functions depend on 
the momenta ^ . The parameter /i remains the renormalization scale while 
the factorization scale is chosen as the hard scale Q. 

The predictive power of the factorized form of the amplitude follows from 
the properties of its individual functions. The jet functions in include all 
coUinear dynamics, and hence all double poles in dimensional regularization, 
which arise from the overlap of infrared and coUinear enhancements in per- 
turbation theory. The matrix of soft functions S^j provides at most a single 
infrared pole per loop and, although it depends on the process kinematics it is 
otherwise completely determined by a set of anomalous dimensions that can 
be computed in the eikonal approximation. The hard-scattering functions /ij^ 
are fully infrared finite. The jet functions and the soft function can be defined 
in terms of specific QCD matrix elements [10,13]. 

The independence of the amplitude Jv[ from the choice of the factorization 
scale Q leads, by the usual connection between factorization and evolution, to 
the renormalization group equation [6] 



din Q 



To shorten the notation, Si,i and hj dependence on pj , will be implicit. 
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where the variations of the jet functions and the soft matrix with the scale Q 
are compensated by variations of the hard function. 

There is stiU considerable freedom in the construction of the jet and soft func- 
tions. For examples, we may shift infrared finite contributions between the jet 
and hard scattering functions, and/or single- logarithmic, process- independent 
soft contributions between the soft matrix and the jets. Specifically, the soft 
matrix Sli is defined only up to any multiple of the identity matrix in the 
color basis space. Below, we will identify a convenient, but by no means unique, 
scheme for the definition of the various functions in (3). 

Of primary importance here is the observation that mixing of color structures 
as a result of soft parton exchange is entirely contained within the matrix F'^l, 
and is therefore summarized by the following solution to the renormalization 
group for the soft function. 



Smf^,a,(/.^),6Upexp 



r,2 



^,as(l^ ),e 



, (6) 



where P stands for path ordering. In the dimensionally-regularized theory, the 
effective coupling should be thought of as expanded in powers of as{fJ.^)- To 
the accuracy we work, this is given by the one-loop form [14] 



a 




(7) 



where 



The integrals arising in the exponential of Eq. (6) are quite simple and for 
comparison to fixed nth-order calculations, we only need to collect all con- 
tributions up to C(Q;^(/i^)). All the one- loop anomalous dimensions in these 
equations have been computed previously in Refs. [10], so that the color mixing 
due to single soft gluon exchange can be predicted (and checked) for specific 
processes. 



In the next section we will use a specific definition of the jet function, based 
on factorization and on the application of Eq. (3) to the simplest color fiow of 
all, the time-like Sudakov form factor. 
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3 Jet and soft functions 



A convenient explicit expression for jet functions may be developed by ap- 
plying the factorization in Eq. (3) to the singlet form factors for quarks and 
gluons, which we denote by To be specific, we consider the quark 

Sudakov form factor. In this case, there is no color mixing at all, and the 
soft anomalous dimension reduces to a number. In the following, we shall ab- 
sorb this number into the evolution of the jets, and reduce the soft function 
to unity by definition. Similarly, given that the anomalous dimension for any 
electro weak vertex vanishes in QCD, we may also reduce the hard function to 
unity in this case. 

We thus define our jet functions by the relation 

(8) 

where is the characteristic momentum transfer, fi is the MS renormalization 
scale (/^^ = nlexp[—e{jE — ln47r)]) and the electromagnetic Sudakov form 
factor in D = 4 — 2e dimensions is given by^ [14,15] 



exp< 



1 fde 



2 J e 

^ 



/CW(«3(/^'),6) 



+ e;l^M-l,aJ^,a,(/i^),e, e 



(9) 



In Eq. (9) the ^ integration is defined order by order in perturbation the- 
ory and as in our discussion of the soft function, the running coupling OLg is 
thought of as being expanded in powers of as{^J^)■ As described in [15,16], the 
functions 7]^, /C'^' and may be read off by comparison to explicit fixed- 
order calculations [17,18]. In particular, /C^ is defined as a series of poles 
in e, while l,Q;s,e) includes the non-singular e-dependence before inte- 

gration. Notice that in Eq. (9) all the Q dependence of the Sudakov form 

^ A similar definition may be given for gluon jets in terms of matrix elements of 
conserved, singlet operators. 



5 



factor is organized by the dimensionally-rcgulatcd running coupling and that 
by construction /C^ = = and 7^ - 7I [15,19]. 

Following the method of Refs. [15] , we verify that 7^ is the familiar Sudakov 
double-logarithmic anomalous dimension known up to two loops [20] to be 



7^ = 2Q(^ 

TT 



(10) 



with Cq = Cf, Cg = Ca and K = Ca (67/18 - ((2)) - Nf{5/9)). 

Before going on to elastic scattering amplitudes, we observe that our choice 
of jet function (8) corresponds to a particular scheme for the factorization of 
the Sudakov amplitude in which the hard-scattering and soft functions are set 
to unity, = 1 in Eq. (5) for the form factor. In this scheme, the hard 

scattering and soft functions of 2 — > 2 amphtudes will be computed with these 
"Sudakov-defined" jet functions. In particular, the matrices Pl^l found in this 
way correspond to those calculated in Refs. [10,13] at one loop. We will return 
to the issue of alternative choices below. 



4 Pole structure for multi-loop hard-scattering amplitudes 



The combination of Eqs. (6) and (8) for the soft and jet functions apphed to 
Eq. (3), allow us to give a useful expression for an arbitrary hard-scattering 
amplitude at two loops and beyond. In particular, the all-orders structure of 
any 2 — > n amplitiude at wide angles is readily understood as the exponenti- 
ation of an appropriate power of the singlet form factors times the expansion 
of the exponentiated soft anomalous dimension. With input to n -|- 1 loops in 
the singlet form factors of quark and gluon, and n loops in the soft anomalous 
dimension, we can predict the full nt/i-next-to-leading poles in dimensional 
regularized amplitudes, in much the same way we predict powers of logarithms 
in threshold or other resummations. 

In particular, on the basis of two- loop calculations for the jet functions and 
of the Sudakov form factor, we can predict almost the entire two-loop single- 
pole term for explicit calculations involving quarks. Conversely, on the basis of 
two-loop calculations we can also determine the soft anomalous dimensions at 
this order, and proceed to a full next-to-next-to-leading pole approximation. 

For ease of comparison, we follow Refs. [9] and denote the vector in color space 
with components A^l in Eq. (2) projected into the orthogonal basis {cj} 

as 
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l-Mf^l)=-ME^({p.},^,«.,e) {cL\^^y (11) 



whose perturbative expansion is 



\M^'^)=Y1 |A^[f(-)l). (12) 

m=0 \ ^ / 



Applying Eqs. (6) - (9) to the factorized form of the amphtude in Eq. (3), 
we easily collect all the singular structure of the one-loop amplitude, using 
jCim ^ 7g^)V(2e), into a function F^'^^^^ as^ 

lM^i^^') = i^'^^^^'(e)|Al[^(°^l) + lM1n^')' (13) 
with the lowest order amplitude 

The function ^[^(^^^(e) in (13) is given in terms of the one-loop coefficients in 
Eqs. (6) and (9) by 



ir[f(i)l(e) = l 



Here, we define 



'tI^ ^)]^^ ^ rim] 



si'^iEe'". I'SAHi^i, (16) 



2 ief 2 ief 



and ^o''"'^''^ = QiQ'^ I = 1, e = 0). The e-poles at O (q;^) are organized by this 
divergent function into those that are color uncorrelated, which are given by 
Sudakov exponentiation through Q and ^k-, and those that are color correlated, 
which are collected by the one-loop soft anomalous dimension matrix pl^^^^J. 

The finite reminder in Eq. (13) is the one- loop hard scattering function of Eq. 
(3), in the notation, 

|AlS») = '.f»(p..f)(c,)„., (17) 



^ In the following, we use the subscript fin to denote finite functions and the sub- 
script ren to denote renormalized functions. 
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At two loops, we can conveniently organize the expansion of Eq. (3) by making 
use of the divergent structure identified at one loop through F^^^^^\ and by 
taking into account that /CK^)] _ -/3o^g^)]/(i6e2) + -f§^^^ / (Ae) and 7^^^' = 



17^'^^. We find 



[f(2)]\ 



where 



with 



4e 



(18) 



(19) 



i,[f(2)](,) 



'[f(2)] 



K^[f(l)] 



1 + r[^(^)i - -rw^^i 




• (20) 



To evaluate the expressions given above, we use 



3 



C 



= C7^|3(^-^C(2) + C(3))C^ 



(21) 



2545 11 , , , A Ca 

+ — C 2 - 13C 3 — 

108 3 ^ ^^7 4 



(22) 



which are found by direct comparison of Eq. (9) with the known fixed-order 
results [15,17,18] for the quark form factor. 



From Eqs. (13) - (15), we see that the bulk of the one-loop singular behavior is 
determined by a combination of the Sudakov amphtudes Al'"^^', from which 
we can find 7^^ and and the one-loop anomalous dimension matrix r'^^^^l. 
Together, these functions control the two-loop singular behavior from down 
to e~^. Information specific to a 2 ^ n process at this order appears first at 
O (e~^) in the contribution of the two-loop anomalous dimension matrix r[^*^^^l, 
which can be determined by comparison to an explicit calculation. We may 
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verify that the expressions above coincide with the proposal of Catani for the 
structure of two-loop hard scattering amplitudes, by direct comparison with 
Ref. [9]. 

As an application, consider the explicit two-loop calculation for the process 
qq qq [2]. To express the amplitude as in Eq. (2), we employ the basis of 
t-channel singlet and octet color exchange (with T" the group generators in 
the fundamental representation) [2,10,13], 

Ci = 5r^rA ^r^TB , C2 = ^ T^^r^ ^t^vb " (23) 

a 

In this basis the one-loop soft anomalous dimension matrix^ is given by 



2CfT 

2{S-U) 



(24) 



N 



where 




r..n^-^j. M^ln^-^j (25, 
and {s, t, u) are the usual Mandelstam variables. 

If we use Eq. (24) and apply the explicit expressions for 7]^*'^'" , 7;^''^'*' and ^0^*'^'" 
given above, wc derive predictions for the poles e^^, e^^ and e^^ that arc in 
complete agreement with the explicit calculation of [2]. We can therefore use 
the term of this calculation to determine the contribution of rt^^^^l for this 
process. More precisely, we can determine its color-uncorrelated part, since 
this is what is given in [2]. 

This direct comparison of the color uncorrelated amplitudes found Refs. [2] for 
qq scattering matrix-elements by diagrammatic evaluation gives the following 
surprisingly simple result 



{vf - |r«) lA^W) = C,(3o (^-^ + ij . (26) 

Thus, comparison of Eqs. (18) - (20) with the explicit calculation shows that 
most of the pole terms at two loops are understandable through exponentia- 

^ It corresponds to the one presented in Eq.(51) of Ref. [10], when we choose the 
scale to be the ji^ = s = {Ia + ^b)^- 
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tion. The remaining coherent part, sensitive to the details of the process, is 
much simplified once exponentiation is taken into account. 

Even this remaining part is subject to further simplification, as is strongly 
suggested by the presence of the overall factor Pq, on the right of Eq. (26). 
In the scheme we have used above, the Sudakov form factor directly deter- 
mines the jet definition. As wc have pointed out, this is equivalent to setting 
^[m^i] _ -[^_ jg easy to show that choosing a nonzero for ^^1**^^(^)1 

produces a corresponding change in Q^^''^^ proportional to /3o . The im- 

portance of the analogous scheme dependence in /cy-resummation has been 
emphasized recently in Ref. [8]. Thus, by using our freedom to modify the 
hard-scattering function in the Sudakov form factor, or equivalently to change 
the normalization of the jet function, we can easily impose the condition 



Given our ability to impose Eq. (27), the uncorrelated part of the single poles 
can be absorbed entirely into appropriately-chosen jet functions. This free- 
dom explains the relations found in Refs. [2,4]. There, on the basis of direct 
calculations it was found that for any 2 — 2 QCD process f with Ug external 
gluons and Ug external quarks and antiquarks 



where ifW^)! is the coefficient of the identity in Iy'^*-^'"(e) on Eq. (20). 

To determine the two-loop anomalous dimension fully, we will need two-loop 
amplitudes as vectors in the basis states, not only in uncorrelated form. These 
have been given so far for gluon-gluon scattering [4], with an expression for 
single-pole terms that is consistent with our results. 



5 Higher orders 

It is clear that the pattern described above extends to higher orders. For 
example, the three-loop hard scattering amplitude has the structure 



\M^m) = FKi)l(e)|A<K2)l) + fK2)1(6)|A<Ki)1) + FK3)](,)|_^[f(o)]^ + |_^[f(3)]^ 




(27) 



(28) 



(29) 
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where 



-[f(3)], 



and 



^[f(i)]/ 



^[f(2)] 

1 r 



^[f(2)], 



ir[f(i)l(e)l'_F[^(2)l(e) 



2F[^(^)](3e)-F[fWl(2e) 



+ ^[^(2)] (3,) _ 1 x[f(2)] (2e) I + ii,[f(3)] (3e) , (30) 



i[f(3)](,) 



Ik I s^'o 
2e2 



[f(3)] ■ 



1 + 



rK3)] 



if! 



(31) 



These expressions enable us to predict the poles in e from sixth order down to 
fourth on the basis of the two-loop Sudakov form factor and the one-loop soft 
anomalous dimension matrices. Beyond this level, the three-loop^ 7^-^^ ^-nd 
two-loop soft anomalous dimension matrices will control all pole structure 
down to double poles, while the color correlations in the single poles will 
determine F^'^)!. 



6 Conclusions and outlook 



We have shown how the factorization properties of hard-scattering amplitudes 
enable us to understand many of the complexities of their infrared poles in 
dimensional regularization. In particular, the universality of exponentiated 
coUinear and soft singularities associated with incoming and outgoing partons 
provides a way of understanding important regularities predicted for and found 
in explicit two-loop calculations. They also predict the structure, and for the 
leading terms the coefficients, of poles to any order. 

The program of calculating two-loop dimensionally-regulated amplitudes is 
a step toward standard model hard-scattering cross sections beyond next-to- 
leading order. In this context, it is useful to note that the factorization in Eq. 
(3) for amplitudes is shared by cross sections in the limit of partonic threshold, 
where the parton center-of-mass energy is just large enough to produce the 
final state. At any fixed order in perturbation theory, leading coUinear and 

^ The fermionic contributions to 7^*^^^' have been calculated recently [21] 
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infrared singularities occur at partonic threshold. Therefore, the exponentiated 
structure of singularities in e for amphtudes will reoccur in the singularities 
of cross sections with radiation. This observation is the basis for threshold 
resummation, which may be a useful guide in organizing the very challenging 
phase space integrals in next-to-next-to-leading order cross sections. 

Here we have explored singularities associated with purely virtual corrections, 
but the simplifications that we have found suggest that a similar analysis may 
be useful for inelastic processes. Leading singularities in cross sections cancel 
between sets of cut diagrams that include the elastic amplitudes at two loops. 
The structure of singularities in these amplitudes, as organized above, may 
help streamline the calculation of cross sections at two loops and beyond. 
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